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Abstract 



Field theoretical renormalization group methods are applied to a simple 
model of a passive scalar quantity advected by the Gaussian non-solenoidal 
("compressible") velocity field with the covariance oc 5{t — t')\x — x'l"^. Con- 
vective range anomalous scaling for the structure functions and various pair 
correlators is established, and the corresponding anomalous exponents are cal- 
culated to the order e"^ of the e expansion. These exponents are non-universal, 
as a result of the degeneracy of the RG fixed point. In contrast to the case of a 
purely solenoidal velocity field (Obukhov-Kraichnan model), the correlation 
functions in the case at hand exhibit nontrivial dependence on both the IR 
and UV characteristic scales, and the anomalous scaling appears already at 
the level of the pair correlator. The powers of the scalar field without deriva- 
tives, whose critical dimensions determine the anomalous exponents, exhibit 
multifractal behavior. The exact solution for the pair correlator is obtained; it 
is in agreement with the result obtained within the £ expansion. The anoma- 
lous exponents for passively advected magnetic fields are also presented in the 
first order of the e expansion. 

PACS number(s): 47.10.+g, 47.25.Cg, 05.40.+j 
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Typeset using REVT^ 



I. INTRODUCTION 



Much attention has been paid recently to a simple model of the passive advection of 
a scalar quantity by a Gaussian short-correlated velocity field, introduced by Obukhov 
and Kraichnan 0, see the papers [3-24] and references therein. The structure functions 
of the scalar field in this model exhibit anomalous scaling behavior, and the corresponding 
anomalous exponents can be calculated explicitly using certain physically motivated "linear 
ansatz" 0, within regular expansions in various small parameters p|- [IO| , p!6| , p2| , and using 
numerical simulations P Jl8|j2T|j23[] . On the other hand, this model provides a good testing 
ground for various concepts and methods of the turbulence theory: closure approximations 
PJ^jTSlI, refined similarity relations [lT2|-p!^, Monte Carlo simulations p!5| , p3| , renormalization 
group |2^, and so on. 

The advection of a passive scalar field 6{x) = ^(t,x) is described by the stochastic 
equation 



dt9 + di{vi9) = uo^9 + f, 



where dt = d/dt, di = d/dxi, vq is the molecular diffusivity coefficient, A is the Laplace 
operator, v(x) is the transverse (owing to the incompressibility) velocity field, and / = f{x) 
is a Gaussian scalar noise with zero mean and correlator 



{f{x)f{x')) = 5{t-t')C{Mr), r^|x-x'|. (1.2) 

The parameter L = M^^ is an integral scale related to the scalar noise, and C{Mr) is some 
function finite as L ^ oo. Without loss of generality, we take C(0) = 1 (the dimensional 
coefficient in ( |1.2|) can be absorbed by appropriate rescaling of the field 9 and noise /). 

In the real problem, the field v(x) satisfies the Navier-Stokes equation. In the simpli- 
fied model considered in [2-8], v(x) obeys a Gaussian distribution with zero average and 
correlator 

{v,{x)v,{x')) = Do J rfkP,,(k) {k' + m')-'/'~'/' exp[ik ■ (x - x')], (1.3) 

where Pjj(k) = 6ij — kikj/k"^ is the transverse projector. A; = |k|, Z^o > is an amplitude 
factor, 1/m is another integral scale, and d is the dimensionality of the x space; < £ < 2 
is a parameter with the real ( "Kolmogorov" ) value e = 4/3. The relations 

Do/uo = go = A' (1.4) 

define the coupling constant ( "charge" ) go and the characteristic ultraviolet (UV) momentum 
scale A. 

The quantities of interest are, in particular, the single-time structure functions 

Sr,{r) = {[9{t,^)-9{ty)r), r=|x-x'|. (1.5) 

In the model ( |l.l[ )-( pr^ , the odd multipoint correlation functions of the scalar field vanish, 
while the even single-time functions satisfy linear partial differential equations [0, see also 
Il5|j7|j2^ . The solution for the pair correlator is obtained explicitly; it shows that the structure 
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function 5*2 is finite for M , m = p[ . The higher-order correlators are not found explicitly, 
but their asymptotic behavior for M — can be extracted from the analysis of the nontrivial 
zero modes of the corresponding differential operators in the limits 1/d ^ e — > 

0J^, or e 2 [^|l6l. It was shown that the structure functions are finite for m = 0, and in 
the convective range A >> 1/r >> M they have the form (up to the notation) 

S2n{r) oc Z}o""r"(2-^) (Mr)^", (1.6) 

with negative anomalous exponents A„, whose first terms of the expansion in 1/d and 
e have the form 

A„ = -2n{n - l)e/{d + 2) + 0{e^) = -2n{n - l)e/d + 0{l/d^). (1.7) 



In the paper |^2|], the field theoretical renormalization group (RG) and operator product 
expansion (OPE) were applied to the model (|1.1| )- (|1.3| ). In the RG approach, the anomalous 



scaling for the structure functions and various pair correlators is established as a consequence 
of the existence in the corresponding operator product expansions of "dangerous" composite 
operators [powers of the local dissipation rate] , whose negative critical dimensions determine 
the anomalous exponents A„. The exponent e plays in the RG approach the role analogous 



to that played by the parameter e = A — d in the RG theory of critical phenomena P5 
The anomalous exponents were calculated in [53] to the order of the e expansion for 



the arbitrary value of and they are in agreement with the first-order results obtained in 
the zero-modes approach The RG approach to the stochastic theory of turbulence is 

reviewed in [EgI. 



In the paper [0, a closure- type approximation for the model (|l.l|)-(p^), the so-called 
linear ansatz, was used to derive simple explicit expression for the anomalous exponents for 
any < £ < 2, (i, and n. Although the predictions of the linear ansatz appear consistent 
with some numerical simulations ||3,23] and exact relations [11,19], they do not agree with 



the results obtained within the zero-modes and RG approaches in the ranges of small e,2 — e 
or 1/d. This disagreement can be related to the fact that these limits have strongly nonlocal 
dynamics in the momentum space, which suggests possible relation between deviations from 
the linear ansatz and locality of the interactions, see discussion in Refs. [p!9| , p3| . [| 

The results of the RG approach are completely reliable and internally consistent for 
small £, but the validity of their extrapolation to the finite values of e is not obvious. Most 
numerical simulations have been limited to two dimensions [§|,^ and have not yet been able 
to cover the small e or large d ranges, in which the reliable analytical results are available. 
Therefore, it is not yet clear whether the anomalous scaling in the small e and finite e 
ranges has the same origin, with the exponents depending continuously on £, or there is a 
"crossover" in the anomalous scaling behavior for some small but finite value of e and these 
ranges should be treated separately. 



^ The small e limit can be treated perturbatively, the effective small parameter equals to the 
reciprocal of the significant range of interactions in the momentum space. This range becomes 
infinite as e goes to zero [p7[] . 
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Another important question is that of the universahty of anomalous exponents. The 
exponents A„ in ( |1.7|) do not depend on the choice of the correlator (|1.2|) and on the specific 
form of the infrared (IR) regularization in the correlator ( |1.3| ). It was argued on phenomeno- 
logical grounds in that the anomalous exponents in the Gaussian model can depend on 
more details of the velocity statistics than the exponent e. The exponents indeed change 
when the function 5{t — t') in the correlator (|1.3|) is replaced by some function with finite 
width, i.e., the velocity has short but finite correlation time |TD|, and when the velocity field 



is taken to be time- independent (see Sec. V of the Ref. [p^] ). 

In this paper, we consider the generalization of the model ( |1.1| )-( ]T?^ ) to the case of a 
non-solenoidal ("compressible") velocity field. In this case, the correlator ( |1.3| ) is replaced 

by 

^ ^ 1^ ^'^:^:;^ explik. (X- (1.8) 

The notation is explained below the Eq. ( |1.3| ); the new quantities are the longitudinal 
projector Qjj(k) = kikj/k"^ and the additional amplitude factor Dq > 0. 



One should not expect that a Gaussian, white- noise model of the type ( |1.1| ), (|1.2D, 
^ provides very good approximation for the real compressible advection; however, it can 
be used to illustrate the important distinctions which exist between the compressible and 
incompressible cases, see e.g. p8|j29[] and references therein. 



The aim of this paper is to give the RG treatment of anomalous scaling with non- 
universal exponents; to compare the results of the e expansion with the nontrivial exact 
exponent, and to present analytic results which probably will be easier to compare with 
numerical simulations than the analogous results for the incompressible case. We apply 
the RG method to the model (|1.1|) , (|1.2| ), ( |1.8|) to establish the existence of the anomalous 
scaling in the convective range and to calculate the corresponding anomalous exponents to 
the second order of the e expansion. We show that the single-time two-point correlation 
functions of the powers of the scalar field in the convective range have the form 

(r(t,x)r(t,x')) oc zy-("+f)/2 A-("+P) (Ar)-^"-^''(Mr)^"+^ r = |x-x'|. (1.9) 

for even n + p and zero otherwise. In addition to e and d, the exponents A„ depend on a 
free parameter, the ratio a = D'q/Dq of the amplitudes in the correlator ( [1.8| ). In the first 
order of the expansion in e they have the form 

/ N an(n — l)de ^, , . 

A. = n -1 + e 2) - J ' + 0{e') 1.10 

2(a — 1 -|- a) 

(the results Ai = — 1 -|- e/2 for any a and A„ = n{—l + e/2) for a = are in fact exact). 
We have also calculated the term of the exponent A„, for any d and a; the result is rather 
cumbersome and is given in Sec. [m| . 

The leading term of the convective range behavior of the structure functions ( |1.5|) in the 
model ( |1.8| ) is completely determined by the contribution (6'^"'); it is obtained from ( p.. 9] ) by 
the substitution n —>■ 2n, p ^ and has the form 

^2„(r) oc z/q-" (M/A)^^". (1.11) 
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It follows from ( |1.9|) that the anomalous scaling in the model ( |1.8|) appears already at 



the level of the pair correlation function. The corresponding exponent A2 is found exactly 
for all < e < 2 from the exact solution for the single-time pair correlator, see Sec. |1[ 

e (a — 1) (d — 1) , , 

Ao = -2 ^ — — (1.12) 

(anomalous scaling for the pair correlator with the exactly known exponent was established 
previously in |^ on the example of a passively advected magnetic field). 

In the language of the RG, the non-universality of the exponents ( p..l(]| ), ( |1.12|) is ex- 



plained by the fact that the fixed point of the RG equations is degenerated: its coordinate 
depends continuously on the ratio a (see Sec. |T|). 

In contradistinction with the model ( p.. 3D , where the anomalous exponents are related 



to the critical dimensions of the composite operators (diOdiO)^ [^, the exponents A„ in 



|1.9|) , (|1 . 1 1|) are determined by the critical dimensions of the monomials the powers of 



the field itself, and these dimensions appear to be nonlinear functions of n, see Sec. 0. 
This explains the difference between the convective range behavior of the model ( p..3|) and 
that of the model ( |1.8| ) and makes the limit Dq rather subtle. 

The model ( [1.8|) remains nontrivial in the case d = 1, where the velocity field becomes 
purely potential. One can hope that the one- dimensional case is more accessible to study 
using numerical simulations, than the lowest-dimensional case d = 2 for (|1.3|) , and it will be 
possible to compare the analytic results (|1.9| )- (|1.12| ) with the numerical estimates (despite 
the fact that the structure functions (|1.11| ) are independent of r, the values of the anomalous 
exponents can be extracted from their dependence on M). In the paper the model ( |1.8|) 



has been studied directly for the one- dimensional case in terms of certain potential function 
for the field 6; the analytic expressions for the anomalous exponents obtained within the 
zero-modes technique have been found to agree with non-perturbative numerical results. 
The relationship between our results and the results of is discussed in Sec. |I|. 



The paper is organized as follows. In Sec. 0, we give the field theoretical formulation 



of the model ( |1.1| ), (|1.2|), ( p,.8| ) and derive exact equations for the response function and 
pair correlator of the scalar field. The explicit solution for the pair correlator is obtained 
and the exact expression (|1.12| ) for the corresponding anomalous exponent is derived. In 



Sec. PH we perform the UV renormalization of the model and derive the corresponding RG 
equations with exactly known RG functions (the /3 function and the anomalous dimension). 
These equations have an IR stable fixed point, which establishes the existence of IR scaling 
with exactly known critical dimensions of the basic fields and parameters of the model. The 



solution of the RG equations for the correlation functions ( |1.9|) is given, which determines 



their dependence on the UV scale. In Sec. 0, the dependence of the correlators on the IR 



scale is studied using the OPE, and the relations (|1.9| ), ( |1.11|) are derived. We also discuss 
briefly the RG approach to the model of passively advected magnetic fields introduced in 
In Sec. 0, we present the calculation of the anomalous exponents in the model ( p..8|) 



to the order of the e expansion. The results obtained are briefly discussed in Sec. W% 
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II. EXACT SOLUTION FOR THE PAIR CORRELATION FUNCTION 



The single-time correlation functions of the field 9 in the models of the type (|1 . 1| ) , ( |1.2D , 
( |1.3| ) or ( |1.8|) satisfy closed linear partial differential equations (see also Refs. PJ^,|^). 
Below we give an alternative derivation of the equation for the pair correlation functions 
based on the field theoretical formulation of the problem. 

The stochastic problem (|1.1|) , ( [L.2|) , (|1.8| ) is equivalent to the field theoretical model of 
the set of three fields $ = {9, 9', v} with action functional 

= 9'De9'/2 + 9' [-dt9 - d(y9) + uqA9] - vD-V/2. (2.1) 



The first four terms in ( [^.1| ) represent the Martin-Siggia-Rose-type action |^-^| for the 
stochastic problem ( |1 . 1|) , (|1.2|) at fixed v, and the last term represents the Gaussian averaging 
over V. Here Dq and are the correlators ( p.. 21 ) and ( |1.8| ), respectively, the required 



integrations over x = {t, x) and summations over the vector indices are understood. 

The formulation (|2.1|) means that statistical averages of random quantities in stochastic 
problem (|1 . 1|) , (|1.2| ), ( |1.8| ) coincide with functional averages with the weight exp5'($), so 
that the generating functionals of total {G{A)) and connected (W^(y4)) Green functions of 
the problem are represented by the functional integral 

G{A) = exp W{A) = j'D^^ exp[S(<l>) + (2.2) 

with arbitrary sources A = A^, A^' , in the linear form 

A<^= f dx[A%x)9{x) + A'^'{x)9'{x) + AJ{x)vi{x)]. 



The model ( |2.1|) corresponds to a standard Feynman diagrammatic technique with the 
triple vertex —9'd{v9) = 9'VjVj9 with vertex factor (in the momentum-frequency represen- 
tation) 

Vj = ikj, (2.3) 

where k is the momentum flowing into the vertex via the field 9'. The bare propagators in 
the momentum-frequency representation have the form 

{99')o = {9'9)l = {~iuj + uoer\ (2.4a) 
= C{k){uj^ + ulk^)-\ (2.4b) 
{9'9')o = 0, (2.4c) 



where C{k) is the Fourier transform of the function C{Mr) from ( |1.2| ) and the bare propa- 
gator (vv)o is given by Eq. (|1.8| ). The parameter go = Dq/uq plays the part of the coupling 
constant in the perturbation theory. The pair correlation functions ($$) of the multicom- 
ponent field $ satisfy standard Dyson equation, which in the component notation reduces 
to the system of two equations, cf. 
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G ^(cJ, k) = —icu + Uok'^ — T.g'g{uJ, k), 



(2.5a) 



D{uj, k) = \G{uj, A;)|2 [C{k) + J:g,e'iuj, k)], (2.5b) 

where G{uj, k) = {66') and D{uj, k) = {66) are the exact response function and pair correlator, 
respectively, and Sg/g, Sg/^/ are self-energy operators represented by the corresponding 1- 
irreducible diagrams; the functions Sgg, Sw in the model ( f^.l] ) vanish identically. 

The feature characteristic of the models like ( p.l|) is that all the skeleton multiloop 
diagrams entering into the self-energy operators Sg/g, Tsqiqi contain effectively closed circuits 
of retarded propagators {66') and therefore vanish (it is also crucial here that the propagator 
(vv)o in ( |1.8D is proportional to the 5 function in time). Therefore the self-energy operators 
in (|2.5|) are given by the single-loop approximation exactly and have the formQ 

y (u^k^- [ [ ^o(fc^-(k-q)Vg^)+i^[,(k-q)Vg^ 

S,,,(.;, k)--J—J ^ + ^2)./2+./2 G{q,u), (2.6a) 



^e'e'{u;,k) = J —J ^ (g2 + ^2)./2+./2 (2.6b) 

where g' = |k — q|. The integrations over u' in the right hand sides of Eqs. (|2.6| ) give 
the single-time response function G{q) = (l/27r) / du' G{q,u') and the single-time pair cor- 
relator D{q) = (l/27r) / du' D{q,u)') (note that the both self-energy operators are in fact 
independent of u). The only contribution to G{q) comes from the bare propagator ( p.4aD , 
which in the t representation is discontinuous at coincident times. Since the correlator ( p..8|) , 
which enters into the single-loop diagram for T^e'e, is symmetric in t and t', the response 
function must be defined at t = t' by half the sum of the limits. This is equivalent to the 
convention G{q) = (l/27r) / du' {—iu' + u^k"^)'^ = 1/2 and gives 

V f M r 1 Do(P-(k-q)Vg^) + D(,(k.q)Vg^ 

^O'siu., k) = (-1/2) J ^ (,2 + ^2)./2+e/2 " (2-7) 

The integration over q in ( p.7|) is performed explicitly: 

Ee^eiiv, k) = -e ^° ~ ]^ ^ ^° J{m), (2.8a) 

where we have written 

u^.- f dq 1 _ Tie/2) m-^ 

^ ' J {27fY {q^ + m^Y/^+^/^ {47i)^/^r{d/2 + e/2)' ^' ' 



^ The single-loop approximation to the Dyson equations in the stirred hydrodynamics is equivalent 
5| to the well-known direct interaction approximation (DIA) . One can say that in the models 



(1.1), (O), (O) or (1.8) the DIA appears to be exact. 
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Equations (|2.5a|) , (|2.8|) give an explicit exact expression for the response function in our 
model; it will be used in Sec. |T| for the exact calculation of the RG functions. Below we 
use the intermediate expression ( |2.7| ). The integration of Eq. ( |2.5b| ) over the frequency uj 
gives a closed equation for the single-time correlator. Using ( p.7|) it can be written in the 
form 

(2.9) 

The function C{k) is supposed to be analytic in k"^, which along with the requirement that 
C(fc = 0) = (so that the Eq. ( p. . 1|) has the form of a conservation law for 6) gives 

C{k) = k^<if{k) (2.10) 

with some function "^{k), or in the coordinate representation C{Mr) = — A\E'(r), where 
\E'(r) vanishes rapidly for r ^ oo. 

In the coordinate representation, Eq. ( |2.9| ) takes the form 

2z/oAD(r) = + Do{6,,A - d,d,){A,D{r)) + D',d,d,{A,D{r)), (2.11) 

where we have written 

A M= f gig,- [exp (iq ■ r) - 1] 



For D'q = 0, Eq. (|2.11|) coincides (up to the notation) with the well-known equation for the 



single-time correlator in the model ( |1.3| ) obtained in 0. 

For < e < 2, the equations ( p.9|) , ( p.l2| ) allow for the limit m — > 0: the possible IR 
divergence of the integrals at q = is suppressed by the vanishing of the expressions in the 
square brackets. In what follows we set m = 0. Then Eq. (|2.12|) gives 



Aij{r) = -Br' 



6,j+er,rj/r' , (2.13a) 



- -r(-e/2) 

{4nY/^2'{d + e)T{d/2 + e/2) ^ ' ^ 



(note that B > 0). Using Eq. ( |2.13| ) and the fact that the function D{r) depends only on 
r = |x — x'l, the differential operators entering into Eq. (|2.11| ) are represented in the form 

didj iAijD{r)) = -B{1 + ^y-'^dr (r^'^-^^/^^+'^dr {r'^'^-^^/^^+'^ D{r))) , (2.14a) 

(5,,A - d,d,){A^D{r)) = {d- iy~% (r''-'+'drD{r)) , (2.14b) 
where dr = d/dr; and for the (i-dimensional Laplace operator one has: 

A^(r) = r^-'^dr (/"^9^^(r)) . (2.14c) 
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It then follows from ( |2.14|) that one integration in Eq. (|2.11| ) is readily performed: one 
can just omit the overall "factor" r^~'^drr'^~^; the integration constant is determined by the 
requirement that the solution have no singularity at the origin (r = 0): 

2PodrD = dr"^ - B{d - l)Dor'drD - B{1 + 6)D'f^r-'^'^-^'>/^^+''>dr (t^'^^+^^/^^+'^d) . (2.15) 

Equation (|2.15|) is rewritten in the form 

dr [(1 + hoiArf)^ D{r)] = [1 + /io(Ar)"]^^^9,^, (2.16) 

where we have denoted 

{d- 1) + a{l + e) 



h. = B 



(2.17a) 



and the exponent ( has the form 



^ = ^/2z/o, 



{d + e) D', 



so that 



C-1 



{d-l)D^ + {l + e) D'o' 

{d - 1) {D', - Do) 
{d-l)Do + {l + e) D'o 



(2.17b) 



(2.18a) 



(2.18b) 



Equation ( |2.16| ) is integrated explicitly; the integration constant is found from the require- 
ment that the solution vanish at infinity (including the special case /iq = 0): 



Dir) 



[1 + ho{kr] 



(2.19) 



For = (so that a = C = 0), the expression (p.l9|) reduces (up to the notation) to the 
well known solution for the purely solenoidal velocity field obtained in p[. Dimensionality 
considerations give ^(r) = M-2^(Mr) with some dimensionless function ip, see ( p.lO|) , so 
that ( p.l9| ) can be rewritten as 



Dir) 



-1 



2z/oM2[l + /io(Ar)=]C jmt 



dy[l + ho{Ky/Myf-^dyi,{y). 



(2.20) 



We are interested in the asymptotic form of the correlator D{r) and the structure function 
^2 oc -D(O) —D{r) in the convective range A >> 1/r >> M, where A is determined by ( p.. 41) . 
From (p:20| ) it then follows 



D{r = 0) ~ C M-^-'^^-^^ A<'^-'\ 



(2.21a) 



where we have used the definitions ( |1.4| ) and ( p.l7|) , and C is completely dimensionless factor 
independent of r, M and A: 
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-hi 



2 JO 

For the correlator D{r) in the region A >> 1/r >> M one obtains: 

D{r) ~ h^^{Kr)-'^D{r = 0). 



(2.21b) 



(2.21c) 



It follows from ( |2.21| ) that D{r) differs from -D(O) by the factor oc (Ar)^'^'' << 1. Therefore, 
the leading contribution to the structure function 5*2 oc D{0) ~ D{r) in the convective range 
is given by the constant term -D(O), while the r dependent contribution determines only a 
vanishing correction. Then the comparison of the expression ( |1.11| ) for n = 1 with the exact 
result (|2.21a|) gives A2 = —2 + e — e(, which along with Eq. (|2.18a 



leads to the exact 



expression ( |1.12|) for the critical dimension A2, announced in the Introduction. 

The expressions ( |2.21|) simplify for d = 1 (and for Dq = D'q and any d), when C = 1? 

(EH): 



see 



Dif 



oc I/n^^M"^ 



(Ar) 



(2.22a) 



D( 



0) oc z/( 



^2.22b) 



The expression ( p.22a| ) agreesQ with the result obtained in |^ directly for d = 1. In 
the language of the papers |5H8|j20[| , the leading non-universal term in (|2.22a|) is related to a 
nontrivial zero mode of the differential operator entering into Eq. ( p. 11 ). We note that the 
anomalous scaling for the pair correlator with the exactly known exponent was established 
previously in [^] on the example of a passively advected magnetic field, as a result of the 
existence of a nontrivial zero mode of the corresponding differential operator. We also note 
that for the purely solenoidal case, the analogous zero mode is independent of r and cancels 
out in the structure function, so that the IR behavior of the latter is determined by the 
universal correction term r^"*^. 

In the subsequent Sections, the asymptotic expressions ( p.21| ) will be generalized to the 
case of higher-order correlators and structure functions. 



III. RENORMALIZATION, RG FUNCTIONS, AND RG EQUATIONS 

The analysis of the UV divergences in a field theoretical model is based on the analysis 
of canonical dimensions. Dynamical models of the type (|2.1| ), in contrast to static models, 
are two-scale, i.e., to each quantity F (a field or a parameter in the action functional) one 
can assign two independent canonical dimensions, the momentum dimension d^p and the 
frequency dimension c?'^, determined from the natural normalization conditions d\ = —d^ = 
1, d'^ = dt^ = 0, d^ = df = 0, d'^ = —df = 1, and from the requirement that each term 



^ In the previous version of this paper, it was erroneously stated that the expression ( 2.22a ) 
disagrees with the result |^^. The authors are thankful to M. Vergassola for pointing this out to 
them. 
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of the action functional be dimensionless (with respect to the momentum and frequency 
dimensions separately), see e.g. i26| , |37|j38[| . Then, based on dp and dp, one can introduce 



the total canonical dimension dp = dp + 2d% (in the free theory, dt oc A). 

The dimensions for the model (p.l|) are given in Table |, including renormalized parame- 
ters, which will be considered later on. From Table | it follows that the model is logarithmic 
(the coupling constant go is dimensionless) at e = 0, and the UV divergences have the form 
of the poles in e in the Green functions. The total dimension dp plays in the theory of 
renormalization of dynamical models the same role as does the conventional (momentum) 
dimension in static problems. The canonical dimensions of an arbitrary 1-irreducible Green 
function F = ($ . . . $)i_ir are given by the relations 

d^ = d-N^d^, (3.1a) 
d'^ = l-N^d%, (3.1b) 

dr = d^ + 2d'^ = d + 2- N^d^, (3.1c) 

where = {Ng, Ngi, N^} are the numbers of corresponding fields entering into the function 
F, and the summation over all types of the fields is implied. The total dimension dr is 
the formal index of the UV divergence. Superficial UV divergences, whose removal requires 
counterterms, can be present only in those functions F for which dr is a non-negative integer. 

Analysis of the divergences should be based on the following auxiliary considerations, 
see 0,^,0,|S|: 

(i) From the explicit form of the vertex and bare propagators in the model ( |2.1|) it follows 
that Ngi — Ng = 2Nq for any 1-irreducible Green function, where A'^o > is the total number 
of the bare propagators {60)q entering into the function (obviously, no diagrams with Nq <Q 
can be constructed). Therefore, the difference Ngi — Ng is an even non- negative integer for 
any nonvanishing function. 

(ii) If for some reason a number of external momenta occur as an overall factor in all the 
diagrams of a given Green function, the real index of divergence d'^ is smaller than d^ by 
the corresponding number of unities (the Green function requires counterterms only if d'^ is 
a non-negative integer). 

In the model (|2.1| ), the derivative d at the vertex 9'd{\9) can be moved onto the field 9' 
using the integration by parts, which decreases the real index of divergence: d'-^ = d^ — Ngi. 
The field 9' enters into the counterterms only in the form of the derivative 89' . 

(iii) A great deal of diagrams in the model ( p.l|) contain effectively closed circuits of 
retarded propagators {99')o and therefore vanish. For example, all the nontrivial diagrams 
of the 1-irreducible function {99'w)i_i^ vanish. 

From the dimensions in Table | we find dr = d + 2 — N^ + Ng — {d + \)Ngi and d'-p = 
{d + 2) (1 — Ng') — A'v + Ng. From these expressions it follows that for any d, superficial 
divergences can only exist in the 1-irreducible functions with Ng/ = 1, N^ = Ng = {dr = 1, 
4 = 0), Ng, = N^ = Ng = 1 {dr = 1, 4 = 0), and Ng, = Ng = 1, N^ = {dr = 2, 
c/'p = 1) (we recall that Ng < Ng,, see (i) above). However, no diagrams can be constructed 
for the first of these functions, while for the second function, all the nontrivial diagrams 



vanish (see (iii) above). As in the case of the purely solenoidal field [^], we are left with 
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the only superficially divergent function {6'6)i^ir; the corresponding counterterm necessarily 
contains the factor of d6' and is therefore reduced to 6'A6. Introduction of this counterterm 
is reproduced by the multiplicative renormalization of the parameters Qq, uq in the action 
functional ( p.lD with the only independent renormalization constant Z^: 

uq = vZy^ (3.2a) 



5(0 = fl-Zi^^g, (3.2b) 

= Z-\ (3.2c) 

Here [i is the renormalization mass in the minimal subtraction scheme (MS), which we 
always use in what follows, g and v are renormalized analogues of the bare parameters g^^ 
and z/q, and Z = Z{g,a,e,d) are the renormalization constants. Their relation in ( p.2c| ) 



results from the absence of renormalization of the contribution with in ( |2.1| ), so that 
Dq = gQi>Q = gfi'^u. No renormalization of the fields and the parameters m, M, a is required, 
i.e., Zq, = 1 for all $ and mo = m, Z^ = 1, etc. The renormalized action functional has the 
form 



Sreni^) = 6' DeO' l2 + e' [-dtO - d{ve) + uZ.Ae] - vD,-V/2, (3.3) 

where the contribution with is expressed in renormalized parameters using (|3.2|) . 

The relation S'($,eo) = Sren{^,e, fi) (where cq is the complete set of bare parameters, 
and e is the set of renormalized parameters) for the generating functional ^^(A) in ( p.2|) 
yields W^(A, cq) = Wreni^, e, /x). We use to denote the differential operation fid^ for fixed 
Co and operate on both sides of this equation with it. This gives the basic RG differential 
equation: 

VRGWren{A,e,fl) =0, (3.4a) 

where I^_rg is the operation expressed in the renormalized variables: 

Vnc =V, + P{g)dg - i,{g)V,, (3.4b) 

where we have written = xdx for any variable x, and the RG functions (the (3 function 
and the anomalous dimension 7) are defined as 

7^(5() = "D^lnZ^, (3.5a) 

(3{g)=V^g = g{-e + ^,). (3.5b) 



The relation between (3 and 7 in ( |3.5b| ) results from the definitions and the relation ( |3.2c| ). 



The renormalization constant Zj, is found from the requirement that the 1-irreducible 
function {6'6)i_ir expressed in renormalized variables be UV finite (i.e., be finite for e ^ 0). 
This requirement determines Z^, up to an UV finite contribution; the latter is fixed by the 
choice of a renormalization scheme. In the MS scheme all renormalization constants have 
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the form "1 + only poles in e." The function G ^ = {9' 6)1^1.,. in our model is known exactly, 
see Eqs. (|2.5a|) , ( p.8|) . Let us substitute (|3.2| ) into Eqs. (|2.5a|) , ( p.8|) and choose to cancel 



the pole in e in the integral J{m). This gives: 

^ d — 1 + a , , 

where we have written Cd = Sd/ {2nY and Sd = /T{d/2) is the surface area of the unit 
sphere in d-dimensional space. Note that the result ( |3.6| ) is exact, i.e., it has no corrections 
of order 5?^, g^, and so on; this is a consequence of the fact that the single-loop approximation 
(^.81) for the response function is exact. Note also that for a = (|3.6D coincides with the 



exact expression for Zy in the "incompressible" case obtained in [p2| . 

For the anomalous dimension '^^{g) = 'D^lnZ^ = (3{g)dghi Z^ from the relations (3.5b) 
and ( |3.6| ) one obtains: 

-eVg \nZy d-l + a 

From ( p.5b|) it then follows that the RG equations of the model have an IR stable fixed point 
[P{9*) = 0, P'{g*) > 0] with the coordinate 

_ 2de 

- Cdid-l + aY ^^-^^ 

The fixed point is degenerated: its coordinate g^, depends continuously on the parameter 
a = D'q/Dq!^ The value of 'juig) at the fixed point is also found exactly: 

l*u = lu{g*) = e. (3.9) 

The solution of the RG equations on the example of the stochastic hydrodynamics is 
discussed in detail in Refs. [^,^; see also for the case of the model (|1.3|); below we 
confine ourselves to the only information we need. 

In general, if some quantity F (a parameter, a field or composite operator) is renormalized 
multiplicatively, F = ZpFren with certain renormalization constant Zp, its critical dimension 



is given by the expression (cf. p6 i p7 i p8i): 

A[F] = A^ = 4 + A^d^ + 7;, (3.10) 

where d'l^ and d% are the corresponding canonical dimensions, 7^ is the value of the anoma- 
lous dimension jrig) = T^n^'^Zp at the fixed point, and = —2 + 7* = —2 + £ is the 
critical dimension of frequency. The critical dimensions of the fields $ in our model are 
found exactly; they are independent of the parameter a and coincide with their analogues 
in the model (O), cf. p2[: 



^ Formally, a can be treated as the second coupling constant. The corresponding beta-function 
f3a = T>^a vanishes identically owing to the fact that a is not renormalized. Therefore, the equation 
/3q = gives no additional constraint on the values of the parameters 5, a at the fixed point. 
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Av = 1 - 



(3.11a) 



Ae = -l+e/2, (3.11b) 



Ae> =d+l-e/2 (3.11c) 



(we recall that the fields in the model ( p.l| ) are not renormalized and therefore 7$ = for 
all $). 

Let G{r) = {Fi{x)F2{x')) be a single-time two-point quantity, for example, the pair 
correlation function of the primary fields $ = {6, 6', v} or some multiplicatively renormaliz- 
able composite operators. The existence of the IR stable fixed point implies that in the IR 
asymptotic region Ar >> 1 and any fixed Mr the function G{r) is found in the form 

G{r) ^ z/q^ A'^«(Ar)-^« ^(Mr), (3.12) 

with certain, as yet unknown, scaling function ^ of the critically dimensionless argument 
Mr. The canonical dimensions (Iq, do and the critical dimension ot the function G{r) 
are equal to the sums of the corresponding dimensions of the quantities Fj. 

Now let us turn to the composite operators of the form 9^{x) entering into the structure 
functions ( |1.5| ) and the correlators ( [1.91 ). 



In general, counterterms to a given operator F are determined by all possible 1- 
irreducible Green functions with one operator F and arbitrary number of primary fields, 
r = {F{x)(^{xi) . . . <^{x2))i-ir- The total canonical dimension (formal index of divergence) 
for such functions is given by 

dr = dp — Nipd^, (3.13) 

with the summation over all types of fields entering into the function. For superficially 
divergent diagrams, the real index d'j^ = d^ — Nqi is a non-negative integer. From Table I 
and Eq. (|3.13|) for the operators O'^i^x) we obtain dp = —n, d^ = —n + Ng — — {d + l)Ng', 



and dp = —n + Ng — — (d + 2)Nqi. From the analysis of the diagrams it follows that the 
total number of the fields 6 entering into the function F can never exceed the number of the 
fields 6 in the operator 6" itself, i.e., Ng < n. Therefore, the divergence can only exist in 
the functions with A'v = 0, Ngi = 0, and arbitrary value of n = Ng, for which dr = d'j^ = 
and the corresponding counterterm has the form It then follows that the operator 6^ is 
renormalized multiplicatively, 6*" = Zn[0"']ren- 

Note an important difference between the case of a purely transversal velocity field ( [L.3|) 
and the general case ( |1.8| ). In the first case, the derivative d at the vertex can be moved 
onto the field 6 owing to the transversality of the velocity field, 9'd{w9) = 9'{wd)9. This 
reduces the real index d'^ by at least one unity, so that d^ becomes strictly negative, see p2 



This means that the operator 6'" requires no counterterms at all, i.e., it is in fact UV finite, 
Zn = l.Q It then follows that the critical dimension of ^"(x) in the model ( |1.3| ) is simply 



^ This "non-renormalization" result can be interpreted as the fact that the scalar field remains a 
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given by the expression ( |3.1CI| ) with no correction from and is therefore reduced to the 
sum of the critical dimensions of the factors ||22||: 



A„ = A[r] = nA[e] = n{-l + e/2). 



(3.14) 



In the general case ( |1.8|) , the constants Z„ are nontrivial, and the simple relation ( p.l4|) is 
no longer valid. The two-loop calculation of the constants Zn is explained in detail in Sec. 
0, and here we only give the two-loop result for the critical dimensions A„ in the model 
O): 



A, 



, , an(n — l)de a(a — l)n(n — l)(d ~ l)e'^ 

n{-i + e/2) - ^ / , + ^ ^; / . . — + 



2{d-l + a) 2{d-l + ay 

a^n{n-l){n-2)dh{d)e^ 
^ 4(d-l + a)2 



where we have denoted 



A;! 



^^'^^~E>^Hd/2 + l)...id/2 + k) 



F(l,l;rf/2 + l;l/4), 



(3.15) 



(3.16) 



and F(. . .) is the hypergeometric series, see . 
In the special case n = 2 one obtains from (|3.15| ): 



Ao = -2 



e{d - l)(a - 1) ^ e^{d-l)aia-l) ^ ^^^3^ 



(d-l + a) 



(d-l + a) 



(3.17) 



Expression ( |3.15 ) is simplified for any integer value of d owing to the fact that the series 
in (|3.16|) reduces then to a finite sum, see [p9| : 



h{d) = 2d 
for any even value of d and 
h{d) = 2d 



d/2 /'_q\A;-2 



(3.18a) 



{d-l)/2 /_q\(d-l)/2-fc 
(_l)(rf-l)/2. 3^/2-2. ^ ^ ^ 

k=l 



2k -1 



(3.18b) 



for any odd value of d, which gives h{d) = 27r/(3v^) for d = 1, h{d) = 41n(4/3) for d = 2, 
and h{d) = 12 — 27rv^ for d = 3. (We note that for ci = 1 and d = 2 the sums in ( p.l8| ) 
contain no terms). The case of a purely potential velocity field is obtained for Dq = const , 



continuous function even in the limit 
ization of the monomials (dOdO) 



or equivalently A — > 00. The nontrivial UV renormal- 



|22| points to the fact that the scalar field is not differentiable, 
i.e., its gradients exist only as distributions. One of the authors (N.V.A) is thankful to G.L.Eyink 
for pointing this out to him, see also Refs. P-H]. 
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Do = or, equivalently, a ^ oo, = g^a = const . From ( p.8|) it then follows that at the 
fixed point = 2de/Cd', the values of the critical dimensions are obtained simply by 
taking the limit a ^ oo in the expressions ( p.l5| ), ( p.l7| ) and have the form 

An = n(-l + 6/2) - n{n - l)de/2 + n{n - l){d - \)e^ 12 + 

- \){n - 2) h[d) de'^/A + 0{e^). (3.19) 

In the special case d = 1 one obtains 

An = —n + ne — n^el2 + n{n — l)(n-2)£V(6v^) + (3.20) 

For the pair correlators of the operators 6'" we obtain from Table I and Eqs. ( p. 12 ), 
&. 

(r(x)r(x')) = z/o"^"+^^/'A-("+^)(Ar)-^"-^''en,p(Mr), (3.21) 

with the dimensions A„ given in ( p.l5| ) and certain scaling functions ^„p(Mr) (for odd n+p 
they vanish). We recall that the representation ( p.21|) holds for Ar >> 1 and any fixed Mr; 
the behavior of the functions ^n,p{Mr) for Mr << 1 (convective range) is studied in the 
subsequent section. 



IV. OPERATOR PRODUCT EXPANSION AND ANOMALOUS SCALING 

The representation ( |3.21| ) for any functions C,n,p{Mr) correspond to IR scaling in the 
region Ar >> 1 and any fixed Mr with definite critical dimensions A„ given in ( 3.15| ). The 



expressions ( |1.9| ) should be understood as certain additional statements about the exphcit 
form of the asymptotic behavior of the functions ^„p(Mr) for Mr —>■ 0. The form of 
the scaling functions ^n,p{Mr) in the representation ( p.21| ) is not determined by the RG 
equations themselves; these functions can be calculated in the form of series in e. However, 
this e expansion is not suitable for the analysis of their behavior for Mr 0, because the 



actual expansion parameter appears to be £ln(Mr) rather than e itself, cf. [p2| , p6| , p8| . In 
contrast to the "large UV logarithms" In(Ar), the summation of these "large IR logarithms" 
is not performed automatically by the solution of the RG equations. 

In the theory of critical phenomena, the asymptotic form of scaling functions for M — 
is studied using the well known Wilson operator product expansion (OPE), see e.g. |]25|; the 
analogue of L = is there the correlation length Vc- This technique is also applied to 
the theory of turbulence, see e.g. P2| , PB| , P5| . 

According to the OPE, the single-time product Fi{x)F2{x') of two renormalized operators 
at X = (x + x')/2 = const , and r = x — x' ^ has the representation 

F,{x)F2{x') = ^C„(r)F,(x,t), (4.1) 

a 

in which the functions Ca are the Wilson coefficients regular in and F^ are all possible 
renormalized local composite operators allowed by symmetry, with definite critical dimen- 
sions Aq. 
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The renormalized correlator {Fi{x)F2{x')) is obtained by averaging ( [4.1|) with the weight 
exp Sr, the quantities (Fa) appear on the right hand side. Their asymptotic behavior for 
M — i> is found from the corresponding RG equations and has the form 

(F„)ocM^". (4.2) 



From the operator product expansion ( [4.1|) we therefore find the following expression for the 
scaling function ^(Mr) in the representation ( |3.12D for the correlator {Fi{x) F2{x')) : 



e(Mr) =5^A„(Mr)^", (4.3) 

a 

with coefficients Aa = Aa{Mr), which are regular in (Mr)^, generated by the Wilson coef- 
ficients in (13). 

We note that for a Galilean invariant product Fi{x)F2{x'), the right hand side of Eq. 
( [4.1| ) can involve any Galilean invariant operator, including tensor operators, whose indices 
are contracted with the analogous indices of the coefficients Ca- Without loss of generality, 
it can be assumed that the expansion is made in irreducible tensors, so that only scalars 
contribute to the correlator {F1F2) because the averages (Fa) for non-scalar irreducible 
tensors vanish. For the same reason, the contributions to the correlator from all operators 
of the form dF with external derivatives vanish owing to translational invariance. 

The leading contributions for Mr are those with the smallest dimension A^, and 
in the e expansions they are those with the smallest da = d[Fa] for e = 0. In the standard 
model (j)^ of the theory of critical behavior one has Aq, = na + 0{e), where no, > is the 
total number of fields and derivatives in F^. The operator F = 1 has the smallest value 
Ha = 0, and it gives a contribution to ([4. 3D which is regular in {Mr^ and has a finite limit 



as Mr 0. The first nontrivial contribution is generated by the operator 0^ with = 2. 
It has the form (Mr)^^*^*^*^-* and only determines a correction, vanishing at Mr 0, to the 
leading term generated by the operator F = 1. 

The distinguishing feature of the models describing turbulence is the existence of "dan- 
gerous" composite operators with negative critical dimensions The contributions 
of the dangerous operators into the operator product expansions lead to a singular behavior 
of the scaling functions on M for Mr -^0. It is obvious from ( p.l5| ) that all the operators 
6"' in the model ( p.l|) are dangerous at least for small e, and the spectrum of their critical 
dimensions is unbounded from below. If all these operators contributed to the OPE like 
( [4.1| ), the analysis of the small M behavior would imply the summation of their contribu- 
tions. Such a summation is indeed required for the case of the different-time correlators 
in the stochastic Navier-Stokes equation, and it establishes the substantial dependence of 
the correlators on M and their superexponential decay as the time differences increase, see 
2^,38]. Fortunately, the problem simplifies for the model ( |2.1| ). 



From the analysis of the diagrams it follows that the number of the fields 9 in the 
operator F^ entering into the right hand sides of the expansions ( |4.1j ) can never exceed the 
total number of the fields 6 in their left hand sides. Therefore, only finite number of operators 
6"' contribute to each operator product expansion, and the asymptotic form of the scaling 
functions is simply determined by the operator 6"' with the lowest critical dimension, i.e., 
with the largest possible number of the fields 6. For the scaling functions C,n,p{Mr) entering 
into the expressions ( |3.21| ) this gives ^„^p(Mr) oc (Mr)^"+p, which lead to the asymptotic 



expression (|1.9|) announced in the Introduction. 
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It is noteworthy that the set of the operators 6^ is "closed with respect to the fusion" 
in the sense that the leading term in the OPE for the pair correlator (6^6^) is given by the 
operator from the same family with the summed index n + m. This fact along with 
the inequality 



A„ + A„ > A 



(4.4) 



which is obvious from the explicit expression ( p.l5|) for small values of e, can be interpreted as 
the statement that the correlations of the scalar field 6 in the model ( ^.11 ) exhibit multifractal 
behavior, see [^0|-^ . In the case of the solenoidal velocity field, the dimension A„ becomes 
linear in n, see (|3.14| ), and the relation ([L9| ) reduces to the so-called "gap scaling" (see [0). 
In this case, the nontrivial multifractal behavior manifests itself in the correlations of the 



dissipation rate rather than in the correlations of the field itself, see [p^j . 

Now let us turn to the structure functions ( |1.5| ) in the convective range Ar >> 1, 
Mr << 1. From the expression (|1.9|) it follows 



2n 



A2„ 



k+p=2n 
k,pj^O 



Cfep(Ar 



I A2n — Afc— Ap 



(4.5) 



where the coefficients Ckp are independent of the scales A, M and the separation r. It is 
obvious from the inequality ( [4.4[ ) that all the contributions in the sum in ( [4.5| ) vanish in the 
region Ar >> 1, so that the leading terms of the structure functions do not depend on r 
and are given by the Eq. ( |1 . 1 1| ) . 

The comparison of the expressions ( |1.9| ) and ( |1.11|) for k = p = 1 with the exact results 
( [^.21cD and ( [^.21a| ) gives A2 = —2 + e — ed^, which along with Eq. ( p^.lSaj ) leads to the exact 
expression (|1.12|) for the critical dimension A2, announced in the Introduction. We note 
that the expression (|3.17|) for A2 obtained within the RG approach is in agreement with the 
corresponding terms of the expansion in e of the exact exponent ( |1.12| ) for all d and a. We 
also note that (|2.21c| ) is consistent with the exact RG result Ag = —1 + e/2, see (|3.11b|) . 

It is seen from ( ^.5| ) that the IR behavior for the structure functions is determined by 
the contributions of the composite operators 6"' to the corresponding OPE. The operators 

obviously do not appear in the naive Taylor expansions of the structure functions ( [L.5| ) 
for r — > 0: the Taylor expansion for the function starts with the monomial {diOdiO)'^. 
However, the operators entering into operator product expansions are not only those which 
appear in the Taylor expansions, but also all possible operators which admix to them in 
renormalization. One can easily check that all the monomials 9"^^ with p < n admix to 
(diOdiO)"' in renormalization. As a result, their contributions appear in the OPE for the 
structure functions and dominate their IR asymptotic behavior. 

The situation changes if the velocity field is purely solenoidal, with the correlator given 
in ( |1.3|) . In this case, the field 6 enters into the vertex in the form of a derivative, 9'd{\9) = 
6'{\'d)6, and therefore only derivatives of 6 can appear in the counterterms to the monomials 
{diddiO)^. Hence, the operators of the form 6'" cannot admix in renormalization to the 
monomials (diOdiO)^ and cannot appear in the OPE for the structure functions ( |1.5| ). This 
means that the contributions of the operators 6^ to the pair correlators ( |1.9| ) cancel out 
in the structure functions, and the IR behavior of the latter is dominated by the operators 
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[diOdiOY] see (|1.6| ). The cancellation becomes possible due to the fact that the dimension A„ 
for a = is a function linear in n, see (|3.14|) , and therefore all the terms in the square brackets 
in ( ^.51 ) are independent of Ar. In this case, the anomalous exponents are determined by 
the critical dimensions of the powers of the operator di6di6; these dimensions are known up 
to the order of the e expansion [E^ . 



For d = 1, the behavior analogous to (|1.6|) in the model ( |1.8|) is demonstrated by the 
structure functions of the field (j){t, x) defined so that 6{t, x) = dx(l){t, x). In this formulation, 
the problem was studied in the paper [^] within the zero-modes approach and using numer- 
ical simulations. The structure functions of the "potential" (j) are not simply related to the 
structure functions of the primary field 6, but they can be derived directly using the RG tech- 
nique. Obviously, the field enters into the vertex in the form of the derivative 6'dx{vdx4>)- 
Therefore, the operators 0"^ are not renormalized and their critical dimensions are given by 
the relations analogous to ( p.l4|) : A[0"] = nA[0], where A[0] = -1 A[6'] = -2 + e/2, see 



Eq. ( |3.11b| ). The structure functions are then given by the expression analogous to 

S2n ^ /^o~"r"(^-^)(Mr)^2", 

where the part of the anomalous exponents is played by the critical dimensions of 
the operators {dx(f)dx4>)"' = 6*^"" given by Eq. ( |3.20| ). In the notation of |]20| we then have 
C2n = n{A — e) + A2n = 2n — ne{2n — 1) — 2Tm{n — l)(2n — 1)£^/3a/3, in agreement with 
the 0{e) result obtained in ||2^ using the zero-modes approach; the exponent C2 = 2 — e is 
exact. 

Let us conclude this section with a brief discussion of the simple model of a passively 
advected magnetic field considered in In this case, both 6 = 9i{x) and the velocity 

are solenoidal vector fields. The velocity field is taken to be Gaussian with the correlator 
( |1.3|) , and the nonlinearity in (|1 . 1|) has the form VjdjOi — OjdjVi. The anomalous scaling in 
this model also appears already for the pair correlator; the corresponding exponent is found 
exactly . 



The RG analysis given above and in is extended directly to this model. It turns out 
that the expressions for the renormalization constant Z^, the RG functions j3 and 7;^, the 
fixed point (7*, and the critical dimension Ag coincide exactly with the expressions (|3.5|)- 



]8|) and ( 3.11b ) for the model ( p.l|) with the substitution a = 0. For the IR asymptotic 



region, the expressions of the form (|3.12| ) are obtained for the correlation functions of various 
composite operators; the corresponding critical dimensions Ap are calculated in the form 
of the e expansions. In particular, for the critical dimensions A2„ of the scalar operators 
6^^" = {OiOiY we have obtained 

A2n = -2n- ^^^;"^^% 0(e^), (4.6) 



and for the special case n = 1 we have 

' d{d-l) 



A2 = -2 -^^^-4^ + 0(6^). (4.7) 



^ In more realistic models of the MHD turbulence the magnetic field indeed behaves as a passive 



vector in the so-called kinetic fixed point of the RG equations, see [p^. 44 1 
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For the dimensions /S.^^ of the second-rank irreducible tensors OiOjO'^"' ^ — SijO'^"' /d we have 

K,^^-2n^ ^^'^'^l-_'^^^-l^;l"-'K o(e^). (4.8) 

The leading terms of the small Mr behavior of the scaling functions are determined by 
the contributions of the scalar operators 6'^'^, and the part of the anomalous exponents is 
played by the dimensions ( [4. 61) . For the special case of the pair correlator it then follows 
{e{x)e{x')) oc {Kr)-'^^o{Mr)^\ In the notation of |30[ we have 7 = A2 - 2A0; from Eqs. 
( pibl) and Q it follows A2 - 2A9 = -e - 2e^{d - 2)/d{d - 1) + 0{e^) for any d, in 
agreement with the exact expression for 7 obtained in 



V. CALCULATION OF THE ANOMALOUS EXPONENTS TO THE ORDER 



In this section we present the two-loop calculation of the critical dimensions A„ of the 
composite operators which determine the anomalous exponents in the expressions ( p..9D , 

The operators 6^ are renormalized multiplicatively, ^" = Zn[0'^]ren (see Sec. |T|). The 
renormalization constants Zn can be found from the requirement that the 1-irreducible 
correlation function 

(r],en(x)^(Xi) . ..e{x2))l-rr = (x)^^(xi ) . . .^(X2))l_„. = Z'^T^ (5.1) 

be UV finite, i.e., have no poles in e, when expressed in renormalized variables using the 
formulas ( |3.2| ). This requirement determines Z^ up to an UV finite part; the choice of the 
finite part depends on the "subtraction scheme" . Most convenient for practical calculations 
is the minimal subtraction (MS) scheme. In the MS scheme, only poles in e are subtracted 
from the divergent expressions, and the renormalization constants have the form "1 + only 
poles in e" . In particular, 

oo oo n 

Z-^ = 1 + E ^k{9)e-^ = 1 + E E ^nke-\ (5.2) 

fc=l n=l k=l 

The coefficients Qnk in our model depend only on the space dimension d and the completely 
dimensionless parameter a; their independence of e is a feature specific for the MS scheme. 
One-loop diagrams generate contributions of order g in ( |5.2| ), two-loop ones generate con- 
tributions of order g"^, and so on. The order of the pole in e does not exceed the number of 
the loops in the diagram. 

The two-loop diagrams of the function r„ required for the calculation of Zn to the order 
g^ and the corresponding symmetry coefficients are given in Table II. The solid lines in the 
diagrams denote the bare propagator {66')q from ( |2.4a| ), the end with a slash correspond to 
the field 9', and the end without a slash correspond to 9; the dashed lines denote the bare 
propagator (|1.8|). Note that the propagator {99)o does not enter into the diagrams for r„. 
The black circle with p > attached "legs" denotes the vertex factor Fp given by the p-fold 
variational derivative Fp = 69^{x)/69{xi) . . . 59{xp). 

Now let us turn to the calculation of the diagrams from Table II. It is sufficient to 
calculate the function r„ in the momentum-frequency representation with all the external 
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momenta and frequencies equal to zero; the IR regularization is then provided by the "mass" 
m from the correlator ( |1.8|) . In what follows, we use the notation 



Rij{k) = DoPijik) + D'oQi^ik) (5.3a) 

and 

5(A;) = (A;2 + m2)-^/2-^/2. (5.3b) 

We also recall the relations Dq = qqUq and a = Dq/Dq. 

The diagram D2 differs from Di only by the insertion of the simplest self-energy diagram 
Sgg/ into one of the two lines {06'). Therefore, the combination Di + 2D2 entering into r„ 
can be easily calculated as a whole: we calculate the single-loop diagram Di with the exact 
propagators {06') instead of the bare propagators {00')q and then expand the result in qq 
to the order Qq. From the exact solution (see Sec. ||) it follows that the propagator {00') 
is obtained from its bare counterpart simply by the replacement 1^0 —>■ i]q, where the exact 
"effective diffusivity" has the form[] 

Vo = J^o + ^ J{m), 

see ( [^.5a| ), (^.8|). Then the "exact" analogue of the diagram Di is given by 

du; f dk S{k)R.,{k)hk, 



27r J (27r)'^ \iuj + r]Qk'^\' 



where we have performed the elementary integration over the frequency and used the isotropy 
of the function S{k). The expansion of the result (|5.4| ) in gives 



go{d — 1 + a)J{m) 



2d 



(5.5) 



The right hand side of Eq. (|5.5| ) is expressed in renormalized variables by the substitution 
9o = 9f^'^Zu^ with the constant Z,^ from (|3.6|), which within our accuracy gives: 



D, + 2D2 = ^^^^ + ^c,/e - /.V(m)] ^ 

= g D^^^ + g"" D^'\ (5.6) 

The diagram is represented by the integral 

_ [ du r duj' f dk f dg^ Rij{q){k + q)i{k + q)jRps{k)kpks ^^j^^^^ . , . 
^''~J 2^J 277 (2^7 (2^ \ico + Mk + q)^\^\ico' + Uok^\^ ^Wb[q), (5.^ 

and the integrations over the frequencies give 



^ See [45 1 for the exact expression for the effective diffusivity in the incompressible case. 
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f dk 


f dq 


4 } 


{27tY J 





a + {1 — a) 



k'^ sin^ {} 
(k + q)2 



S{k)S{q), 



(5.8) 



where "i? is the angle between the vectors k and q, so that k • q = A;g cos 'd. The symmetry of 
the integral ( p.8| ) in k and q allows one to perform the substitution k"^ — * (k + q)^/2 — k • q 
in the integrand, which gives: 



D3 = -^J^{m) + ^ ^ ^ - 2J2M], 



4 



where we have written 



dk 



"^"^ sin2^9 5(fc)5(g) 



and 



Jo Cm) 



t/k 



(27r)'^y (27r)'^ 

(iq k ■ q sin^ 'd 



S{k)S{q). 



{2'kYJ {2'kY (k + q)2 
The integral in ( ^.101) can be easily expressed via J{m): 

roc roo r 

Ji{m) = Cj / dkk'^-^ / dqq'^-^ / dn sin^ ^ S{k)S{q) 
Jo Jo J 



(5.9) 



(5.10) 



(5.11) 



j'^(m) I dnsin^ i9 = — - — j'^{m) 



(5.12) 



with the coefficient Cd from ( p.6|) . Here and below / cin denotes the integral over the d- 
dimensional sphere, normalized with respect to its area, so that / dn 1 = 1 and / dwsiv? 'd = 
{d — l)/d. For the integral (|5.11|) one has: 



J2(m) 



dk / dq 



/c'^g"' cos -(9 sin^ ^9 , 
dn— — S {k)S[q) 



2Cj 



dk dq dn 



k"^ + q^ + 2kq cost? 



+ g2 + 2kqcos'd 



Sik)Siq), 



(5.13) 



where we have used the symmetry of the integrand and integration area in k and q. 

In order to find the renormalization constant, we need not the entire exact expression 
( p.l3| ) for the integral J2(m), we rather need its UV divergent part. The simple power 
counting shows that the UV divergence of the integral ( |5.13 ) is generated by the region 
in which the both integration momenta k and q are large. Therefore, the integral ( [5.13| ) 
contains only a first-order pole in e, and the coefficient in 1/e does not change when the 
integration area [0, 00] for the momentum k is restricted from below by some finite limit, for 
example, [m, 00]. Furthermore, the IR regularization of the integral is then provided by this 
finite lower limit, and one can simply set m = in the functions S{k), S{q), which gives: 



J2(m) 



dk / dq 
Jo 



dn 



k ^q cos sin^ 'd 
_j_ g2 _|_ 2kq cos "i? 



(5.14) 
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Here and below ~ means the equality up to the terms finite for e — 0. From the dimension- 
ality considerations, it is obvious that J2{fn) = where f{e) contains a first-order 
pole in e. It then follows that 



(5.15) 



(we recall the notation T>m = md/dm). The representation ( ^.15| ) allows one to get rid of 
the integration over k in (|5l^) : 



J2(m) ^ Cj 



m 



'2e „i 



e Jo 



dx I dn 



X ^ cos?? sin 
1 + + 2x cost? ' 



(5.16) 



where we have performed the substitution q = mx. The pole in ( |5.1(i| ) is isolated explicitly, 
the integral is UV convergent and one can set e = in the integrand: 



JM ~ Ci 



m 



-2e „i 



e Jo 



dx / dn 



cos {} sin^ {} 
1 + -|- 2x cost? 



(5.17) 



The integrations in (|5.171 ) are performed explicitly: 

-2e 



Jn(m) ~ C: 



m 



2e 



du t? cos 'd sin t? = 



m-^'{l-d) 



2ec/2 



(5.18) 



Combining the expressions (|5.9| ), ( |5.12| ), and ( |5.18| ) we obtain: 



D, = J\m)gl 



a(l — a){d — 1) 
8d 



a{l — a){d ~ 1) m' 
8ed^ 



-2e 



(5.19) 



Within our accuracy, the renormalization of the expression ( |5.19| ) is reduced to the substi- 
tution go gn'^, which gives: 



^3 = J\rn)^l 



2e 



a^g"^ ag'^{l — a){d — I) 



2^2 «(1 -a)(rf- l)(/i/m)2= 



(5.20) 



Now let us turn to the diagram It is given by the expression 



duj f duj' 



dk 



dq 



27T J 271 J {2nY J {27iY 
Rij(k)ki{k + q)jRps{k)qpqsS{k)S{q) _ 
{iuo + iyok^)\iuj' + z/og2|2(_i(cu + 00') + z/o(k + q)2) ~ 

= ^[j2(m) + J3(m)], 



(5.21) 



where we have performed the integrations over the frequencies and made use of the symmetry 
in k and q; the integral J^ijn) is given by: 



dk f dq k-qS{k)S{q) 



{2ttY J {27tY k^ + q^ + k-q' 



(5.22) 
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Proceeding as for the integral J2{fn) above, we arrive at the expression 

-2e ^1 ^ cost? 



,L(m) ~ C: 



m 



e Jo 



dx / dn 



1 + + X cost? ' 



(5.23) 



which is analogous to the expression ( |5.17 ) for J2{m). In contrast to ( 5.17 ), after the inte- 
gration over X in ( |5.23| ) we arrive at the integral over the angles which cannot be calculated 
explicitly. We rather expand the integrand in (|5.23| ) in cos?9: 







dx / dn 



cos{} 



1 + + X cos 



1 f cos^ ^ 
dx I dn y ^ 

k=0 



1 + x^ 



X COS?? 
" 1 + X2 



(5.24) 



and use the formulas 



dn cos^*^ i!} 



{2k -1] 



dx 



d{d + 2) . . .{d + 2k -2)' 
J dn cos^'^+i ^ = 0, 



;i + x2)2fc+2 4(2A; + 1)!' 
For the series in ( [5. 241 ) this gives (we omit an overall minus sign): 
(2A;+ 1)!!(A;!)2 1 ^ k\ 



(5.25) 



-j^ oo 



oo 

77j 5Z 



M ^ {2k +l)\{d + 2)...{d + 2k) Ad A^{d/2 + 1) . . . {d/2 + k) 



h{d)/Ad, (5.26) 



with the function h{d) entering into the expressions (|3.15|) - (p.l8 



Combining the expressions ( p.21| ), ( ^.23| ), ( |5.26| ), and performing the replacement qq 
gH^ we obtain 



n 72/ ^ 2.«^^^ a^g^h{d)Cl{ii/mf' 
Di = J [m)jj — 



The diagram is simply given by 



32de 



2^2 



(5.27) 



(5.28) 



and Dq contains effectively a closed circuit of retarded propagators and vanishes identically. 
Therefore, the function r„ in the two-loop order of the renormalized perturbation theory 
has the form 

r„ = 1 + "^^{D, + 2D2 + D,) + n{n - l)(n - 2)D, + <n - l){n - 2){n - 3) 



8 



(5.29) 



with the symmetry coefficients from Table II and the explicit expressions for Di given in 
(PD, (iH), (113), and {Wm . 
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Within our accuracy, the renormahzation constant (|5.2|) has the form 



^„ =1 + — + + + 

and the requirement that the function ( [5 .11 ) be UV finite in the first order in g gives: 

o^\\9 , nin - 1) 



(5.30) 



+ 



= UV finite, 



(5.31) 



with the coefficient D'^^^ defined in ( p.6|) . The expansion in e of the integral J(m) from ( p.8| 
entering into the expressions for Di has the form 



fi J[m) = — 



/^(l)-^(d/2) , , , , 



+ 0(e) 



(5.32) 



where '^/'(-z) = d\nT{z)/dz. From ( ^.6| ), ( p.31| ), and the first term of the expansion ( ^.32| ) 
one obtains: 



Oil = —an{n — l)Cd/4. 
The UV finiteness of the function (|5.1|) in the order g"^ imphes: 



(5.33) 



0215^^ , «225'^ , aii9 n{n - 1) 



nin — 11 



e e 2 2 



which along with the expressions (|5^), ( p. 201 ), ( |5.27| ), ( p. 2^ ), ( |5.33| ) and the expansion (|5.32|) 
yields: 



«2l/C^ 



n{n — l)a(a — l)(c? — 1) n{n — \)(n — 2)a'^ h{d) 



+ 



32d 



(5.35a) 



022 /Cj 



a n {n — ly a^n{n — 1) n{n — l)a{a — l){d — 1) 

16 8 ^ Ted 

a^n{n-l){n-2) a^n{n - l){n - 2){n - 3) 



a n{n — 1) 
32 



32 

[an{n~ 1) - 2{a + d - l)/d]. 



(5.35b) 



We note that the 0(1) terms of the expansion ( |5.32| ) cancel out in the Eq. ( p.34| ) and 
therefore give no contribution to the coefficients Ojj. 

For the corresponding anomalous dimension 7„ = In Z„ we have: 

7„ = In Zn = (3{g)dg In Z„, = [s + iM]'^^ In ^n, (5.36) 
with the RG functions I3{g) and '^u{g) from (3^). Within our accuracy (|36D yields 
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(5.37) 



and using the explicit expressions (|5.33|) , ( |5.35| ) one obtains: 

—an{n — l)u n{n — l)a{a — l){d — l)u^ n{n — l){n — 2)a'^h{d)u'^ 



In 



+ 



(5.38) 



where u = gC^- It follows from the explicit expressions (|3.7|), ( |5.33| ), and ( p. 35] ) that the 
coefficient in 1/e in the expression (|5.371 ) for 7„ vanishes: —ang'y^ + 2a22g'^ — (aii5')^ = 0. 
This is a manifestation of the general fact that the function 7^ is UV finite, i.e., it has no 
poles in e. 

Substituting the anomalous dimension ( |5.38| ) into the expression (3.1(J) and performing 
the replacement g ^ g* with g^, from ( p.8| ), we arrive at the desired expression ( p.l5| ) for the 
critical dimension of the composite operator 6'". 

It is worth noting that the case c? = 1 is exceptional in the sense that "there are no angles 
in one dimension." We have performed all the calculation directly in c? = 1 and checked that 
the one-dimensional exponents are indeed obtained from the general expressions like ( p. 19 ) 
by the substitution d = 1. 



VI. DISCUSSION AND CONCLUSION 

We have apphed the RG and OPE methods to the simple model (|1 . 1|) , (|1.2D, (|I.8D, which 
describes the advection of a passive scalar by the non-solenoidal ("compressible") velocity 
field, decorrelated in time and self-similar in space. We have shown that the correlation 
functions of the scalar field in the convective range exhibit anomalous scaling behavior; the 
corresponding anomalous exponents have been calculated to the second order of the e ex- 
pansion (the two- loop approximation), see ( p.l5| )-( [3.2CI| ). They depend on a free parameter. 



the ratio a = D'q/Dq of the amplitudes in the transversal and longitudinal parts of the 
velocity correlator, and are in this sense non-universal. In the language of the RG, the non- 
universality of the exponents is related to the fact that the fixed point of the RG equations 
is degenerated: its coordinate depends continuously on a. 

In contrast to the model ( |1.3| ), where the anomalous exponents are determined by the 
critical dimensions of the composite operators {di6di6)"', the exponents in the model ( p..8| ) 
are related to the critical dimensions of the monomials 6'", the powers of the field itself, 
and these dimensions appear to be nonlinear functions of n. This explains the important 
difference between the anomalous scaling behavior of the model ( |1.3|) and that of the model 
( |1.8|) : in the latter, the correlation functions in the convective range depend substantially on 
both the IR and UV characteristic scales, and the structure functions are independent of the 
separation r = |x — x'|. The monomials 6^ in the model ( |1.8| ) also provide an example of the 
power field operators without derivatives, whose correlation functions exhibit multifractal 
behavior (another interesting example is the field theoretical model of a growth process 
considered in |^2[). Analogous behavior is demonstrated by the model of a magnetic field, 
advected passively by the incompressible Gaussian velocity; the corresponding anomalous 
exponents are calculated to the order e (e^ for the pair correlator). 
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The anomalous exponent for the pair correlation function has been found exactly for all 
< e < 2. Its expansion in e coincides with the result obtained using the RG for all values 
of the space dimensionality d and ratio a. The agreement between the exact exponent for 
the pair correlation function and the first two terms of the corresponding e expansion is also 
estabhshed for a passively advected magnetic field. 

These facts support strongly the applicability of the RG technique and the e expansion to 
the problem of anomalous scaling for the finite values of e, at least for low-order correlation 
functions. 

We note that the series in e for all known exact exponents in the rapid-change models 
have finite radii of convergence, a rare thing for field theoretical models. In the language 
of the field theory, this is related to the fact that in the rapid-change models, there is no 
factorial growth of the number of diagrams in higher orders of the perturbation theory (a 
great deal of diagrams indeed vanish owing to retardation, see discussion in Sec. H). In its 
turn, this fact suggests that the series in e for the unknown exponents (for example, the 
anomalous exponents in the original Obukhov-Kraichnan model) can also be convergent. 



It should also be noted that the asymptotic expressions (|1.9|), (|1.11|) result from the 



fact that the critical dimensions A„ are negative, and that the modulus |A„| increases 
monotonously with n. This is obviously so within the e expansion, in which the sign and 
the n dependence of the dimensions are determined by the first-order terms ( 1.10| ), ( [4. 61) , 
while the higher-order terms are treated as small corrections. However, for finite values of 
e the higher-order terms can, in principle, change these features of the dimensions. Indeed, 
the contribution in the second-order approximation for A„ is positive, see e.g. ( p. 151) , 
so that A„ also becomes positive, provided n is large enough. Of course, this conclusion is 
based on the second-order approximation of the e expansion and is therefore not definitive: 
the higher-order terms of the e expansion contain additional powers of n, so that the actual 
expansion parameter appears to be en rather than e itself, cf. p2|j4^ . Therefore, the correct 
analysis of the large n behavior of the anomalous exponents requires resummation of the e 
expansions with the additional condition that en ~ 1. This is clearly not a simple problem 
and it requires considerable improvement of the existing technique. 
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TABLES 

TABLE L Canonical dimensions of the fields and parameters in the model (2.f). 



F 


e 


e' 


V 






m, M, 11, A 


90 


g,a 


4 





d 


-1 


-2 


-2+£ 


1 


£ 







-1/2 


1/2 


1 


1 


1 











dp 


-1 


d+1 


1 





£ 


1 


£ 






TABLE II. The diagrams of the 1-irreducible Green function {9'^{x)9{xi) . . .6{x2)) in the 
two-loop approximation. 

Diagram Symmetry coefHcient 

n(n - l)/2 

n(n — 1) 
n(n - l)/2 
n{n-l){n-2) 
n(n-l)(n-2)(n-3)/8 
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